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' We develop a model to describe odd parity baryon resonances generated dynamically through a 

^SJ , unitary baryon-meson coupled-channels approach. The scheme applies to channels with light and/or 

— ^ ■ heavy quark content. Distinct features of the model are that, i) the interaction is an S-wave contact 

' one, ii) it reduces to the SU(3) Weinberg- Tomozawa Hamiltonian when light pseudoscalar mesons are 



involved, thus, respecting chiral symmetry, iii) spin-flavor in preserved in the light quark sector, and 
iv) heavy quark spin symmetry is fulfilled in the heavy quark sector. In particular, baryon-meson 
C~ ' states with different content in c or in c do not mix. The model is a minimal one and it contains 

no free parameters. In this work, we focus on baryon resonances with hidden-charm (at least one c 
and one c quarks). We analyze several possible sectors and, for the sector with zero net charm, we 
'(—( '■ write down the most general Lagrangian consistent with SU(3) and heavy quark spin symmetry. We 

^ ^ explicitly study the A'^ and A states, which are produced from the 5-wave interaction of pseudoscalar 

I , and vector mesons with 1/2+ and 3/2+ baryons within the charmless and strangeless hidden charm 

sector. We predict seven odd parity A'^-like and five A-Iike states with masses around 4 GeV, most of 
them as bound states. These states form heavy-quark spin multiplets, which are almost degenerate 
in mass. The predicted new resonances definitely cannot be accommodated by quark models with 
three constituent quarks and they might be looked for in the forthcoming PANDA experiment at 
the future FAIR facility. 
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I. INTRODUCTION 

The possible observation of new states with the charm degree of freedom has attracted a lot of attention over the 
past years in connection with many experiments such as CLEO, Belle, BaBar, and others [1-36]. Moreover, the future 
PANDA and CBM experiments at FAIR facility of GSI [37, 38] will aim at obtaining data in the heavy-flavor sector, 
thus opening the possibility for observation of new exotic states with quantum numbers of charm and strangeness. 
In this regard, a clear goal would be to understand the nature of these states, and in particular, whether they can 
be described with the usual three-quark baryon or quark-antiquark meson interpretation, or qualify better as hadron 
molecules within a baryon-meson coupled-channels description. 

Unitarized coupled-channels approaches have shown to be very successful in describing some of the existing exper- 
imental data. These schemes include, for example, models based on the chiral perturbation amplitudes for S-wave 
scattering of 0~ octet Goldstone bosons off baryons of the nucleon 1/2+ multiplet [39-60]. Recently the charm 
degree of freedom has been incorporated in these models and several experimental states have been described as 
dynamically-generated baryon molecules [61-77]. This is the case, for example, of the Ac(2595), which is the charm 
sector counterpart of the A(1405). Some of these approaches are based on a bare baryon-meson interaction saturated 
with the t-channel exchange of vector mesons between pseudoscalar mesons and baryons [61-70], others make use of 
the Jiilich meson-exchange model [71-73] or some rely on the hidden gauge formalism [74-77]. 

Nevertheless, these models do not explicitly incorporate heavy-quark spin symmetry (HQSS) [78-80], and therefore, 
it is unclear whether this symmetry is respected. HQSS is a QCD symmetry that appears when the quark masses, such 
as the charm mass, become larger than the typical confinement scale. HQSS predicts that all types of spin interactions 
involving heavy quarks vanish for infinitely massive quarks. Thus, HQSS connects vector and pseudoscalar mesons 
containing charmed quarks. On the other hand, chiral symmetry fixes the lowest order interaction between Goldstone 
bosons and other hadrons in a model independent way; this is the Weinberg- Tomozawa (WT) interaction. Thus, it is 
appealing to have a predictive model for four flavors including all basic hadrons (pseudoscalar and vector mesons, and 
1/2+ and 3/2+ baryons) which reduces to the WT interaction in the sector where Goldstone bosons are involved and 
which incorporates HQSS in the sector where charm quarks participate. This model was developed in Ref. [81-83], 
following the steps of the SU(6) approach in the light sector of Refs. [84-87]. In these works, several resonances have 
been analyzed and compared to experimental states, such as the S'-wave states with charm C = 1, 2, 3 [81, 83] together 
with C = —1 states [82]. Also this scheme has been recently extended to incorporate the bottom degree of freedom 
[88] in order to study the nature of the newly discovered Af,(5912) and AJ(5920) resonances [89] as possible molecular 
states. 

In this paper we aim at continuing those studies on dynamically-generated baryon resonances using HQSS con- 
straints. We will discuss extensively the details of the model and how heavy-quark spin symmetry is implemented. 
The model respects spin-flavor symmetry in the light sector and HQSS in the heavy sector, and it reduces to SU(3) 
WT in the light sector respecting chiral symmetry. Moreover, we will focus on the dynamical generation of hidden 
charmed states. The coupled-channels in the hidden charm sectors are characterized by containing equal number of 
c- and c- quarks. As we shall discuss, HQSS does not mix sectors with different number of c- or c- quarks. Thus, 
this model has a symmetry SU(6) x HQSS, with HQSS = SUc(2) x SUe(2) x Uc(l) x Ue(l). We wiU pay a special 
attention to the charmless (C = 0) and strangeless {S = 0) sector. Recent works [74-77] predict the existence of a few 
nuclcon-like states with masses around 4 GcV which result from the baryon-meson scattering in this hidden charm 
sector. In this paper we will analyze these results within our model, and predict the existence of several odd parity 
A- and A'^- like bound states with various spins. These resonances can be organized in heavy-quark spin multiplets, 
which are almost degenerate in mass and that can be subject to experimental detection. 

The predicted new resonances might be subject to experimental detection in the forthcoming PANDA/FAIR ex- 
periment. If conflrmed, they definitely cannot be accommodated by quark models with three constituent quarks. 

The paper is organized as follows. In Sec. II we present the WT interaction implementing heavy-quark spin 
symmetry, and analyze the different hidden charm sectors, classified according to their charmed content. In Sec. Ill 
we introduce the unitarized coupled-channels approach used throughout this work. Sec. IV is devoted to present 
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our results and in Sec. V we summarize the conclusions. In Appendix A we give details for the construction of the 
meson and baryon tensors and the computation of the different matrix elements of the interaction. The tables of the 
interaction matrices for the different baryon-meson channels are collected in Appendix B. 



II. THE MODEL 



A. Weinberg- Tomozawa interaction 

The theoretical model wc use has been developed in previous works for baryon-meson sectors involving light and/or 
heavy quarks, but not for those with hidden charm, so we devote this section to fully specify the model. 

The guiding principle is to blend several well established hadronic symmetries in a model as simple as possible. 
Specifically, to comply with chiral symmetry, S\J{Nf)l x S\]{Np)ji for Np flavors, we require the interaction to 
reproduce the Weinberg- Tomozawa ( WT) Hamiltonian [90-92] , a contact 5- wave interaction, when light pscudoscalar 
mesons are involved. The low energy interaction of soft pseudo-Nambu-Goldstone bosons of the spontaneously broken 
chiral symmetry off any (flavored) target takes the WT universal form^ 

^WT = ^2 J>J^, i^l,...,N]r-l, (2.1) 

where / is the decay constant of the pseudoscalar meson (~ 93MeV), and Jp, are the S\J{Nf) group generators 
(with the standard normalization fiufjki = NpSij, where fijk are the structure constants) for pseudoscalar meson 
and target, respectively. Further, K{s) = fc" + k"^ represents the sum of the incoming and outgoing energies of the 
meson. In the center-of mass frame (CM) 

where ^/s is the total CM energy, M the mass of the target, and m the mass of the pseudo-Nambu-Goldstone meson. 

VwT is the tree level on-shell interaction. The normalization we use is such that the corresponding T-matrix for 
elastic scattering is related to the scattering amplitude by 

1 IM p"^^^ — 1 



where k is the CM momentum.^ 
For three flavors and baryons 
form (we assume exact SU(3) symmetry for simplicity) [48] 



For three flavors and baryons in the (nucleon) octet, the Hamiltonian density of the WT interaction takes the 



Hwt(x) = : tr(B7'^[[^,a^0],B]) : [Np = 3), (2.4) 



where cj) and B are the meson and baryon matrices in the adjoint representation of SU(3). On account of the trace 
cyclic property, and neglecting the meson momentum, the WT interaction can be recast as 

^■.tr{[^,d^4']{B\B}): (2.5) 
From the symmetry group point of view, the interaction in Eq. (2.5) is, schematically, 

^WT = ((M^ ® M)adjoint,antisymmotric ® (^^ ® ^)adjoint) ^i„gl^j , (2.6) 

where adjoint and singlet refers to the adjoint and singlet representations of the flavor group. This corresponds with 
the structure in Eq. (2.1). Af, M\ B and B^ refer to the matrices of annihilation and creation operators of meson 
and baryons (see below). 



^ An extra factor 1/2 is to be added if the projectile and target are identical particles. 

2 T = iM = — T, M and T being the amplitudes defined in [93] and [94] , respectively. 2MT equals M in [95] . 
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It is worth noticing that the eigenvalues of the relevant operator in Eq. (2.1), 2JpJ^, can be written using the 
quadratic Casimir operator of S\J{Np) (see e.g. [96]) 



where ^, fip and fir are, respectively, the SU(A^i?) irreducible representations (irreps) of the system, the Nanibu- 
Goldstone boson (i.e., the adjoint representation) and the target. We use the normalization C2 (adjoint; TV) = N for 
SU(Af). In our convention a positive eigenvalue indicates repulsion and a negative one attraction. 



Next we turn to spin-flavor (SF) symmetry, SU(2iVF), [97-99]. This symmetry has been phenomenologically 
successful in the classification of lowest-lying hadrons as well as in uncovering regularities present in the masses and 
other hadron properties [100, 101]. This is particularly true for baryons, a fact that can be understood from the large 
Nc (number of colors) limit of QCD. In that limit SF becomes exact for the baryon sector [102]. As for mesons, the 
lowest-lying states can also be classified quite naturally according to SF multiplets, but the symmetry works worse 
for the meson spectrum. A prime example of this is provided by the pion and rho mesons. They belong to the same 
multiplet of SU(6) and this would require these two mesons to be approximately degenerated in mass. Also, the pion 
is a collective state identified as the pseudo-Nambu-Goldstone boson of the spontaneously broken approximate chiral 
symmetry, whereas the rho meson mass fits well with a pair of constituent quark-antiquark. Vector dominance also 
suggests that the rho meson should belong to the same chiral representation of the vector current (8, 1) -I- (1, 8), which 
is different from the chiral representation of the pion, (3,3*) -t- (3*, 3). The apparent conflict was solved by Caldi 
and Pagels in two insightful papers [103, 104], where a number of related puzzles are clarified. These authors noted 
that chiral and SF symmetries are compatible, as they can be regarded as subgroups of a larger symmetry group, 
'S{]{2Nf)l X SU(2A^i?)ii;, actually a realization of the Feynman-Gell-Mann-Zweig algebra [105]. In their solution, 
the SF extended chiral symmetry is spontaneously broken, the rho meson being a dormant Goldstone boson of this 
breaking. The collective nature of the rho meson has been confirmed in lattice QCD [106]. As it is well known, exact 
SF invariance is not compatible with exact relativistic invariance [107]. In the Caldi-Pagels scenario, vector mesons 
acquire mass through SF-breaking relativistic corrections which restore Poincare invariance. 

While not an immediate consequence of the QCD Lagrangian, SF symmetry emerges in some limits, such as large 
Nc for baryons, as already noted, and partially also in heavy quark limits. We will turn to this point below. The lack 
of exact relativistic invariance is not unusual in other treatments involving hadron or quark interactions to form new 
hadrons, either bound states or resonances. A good example would be the successful "relativized quark model" of 
Isgur and coworkers [108, 109]. Similarly, in our approach the breaking of relativistic invariance is very mild: the spin 
is treated as just an internal label (as another kind of flavor) and our fields are effectively spinless as regards to the 
kinematics. So we have fully relativistic kinematics with SF as a purely internal symmetry. This entails the following. 
Because angular momentum can be transferred between orbital and spin components, and spin is not conserved under 
Lorentz transformations, a strict relativistic treatment requires fields with different spin to behave differently. In turn 
this yields differences in the off-shell propagators for each spin and this breaks strict SF symmetry. To the extent 
that we consider on-shell particles the various fields behave in the same way, and only the off-shell baryon-meson 
propagator (or loop function) would depend on the spin of the particles involved. We disregard this effect.^ On the 
other hand, because we consider pure S-wave interactions, the spin is separately conserved, and also, near threshold, 
the case of interest to us, the specifically relativistic properties of the spin become irrelevant. It should be noted 
that even approaches with formally relativistic Lagrangians are in practice subject to simplifying approximations in 
vertices and propagators which break exact relativistic invariance without relevant phenomenological implications. 
Also we remark that SF invariance is just our starting point for modeling the interaction. Modifications will be 
introduced below to account for other established properties of QCD, and more importantly, we use physical values 
for hadron masses and meson decay constants in our kinematics. 

The compatibility between SF and chiral symmetries implies that the WT interaction can be extended to enjoy SF 
invariance, and this can be done in a unique way [84]. For the on-shell vertex the extension is simply 



(2J|,J^)^ = C2(Ai; A^f) - C2{iSP-.Nf) - C2{flT-:NF), 



(2.7) 



B. Spin- flavor extended Weinberg- Tomozawa interaction 



WT — 



K{s) 



z = 1, . . . , (2A^f )2 - 1, 



(2.8) 



4/2 



However, any such spin dependence would not be easy to extract from phenomenology as it will bo masked by the intrinsic ambiguity 
of the loop function, which has to be renormalized using some phenomenological proscription. 
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where Jlj and are the S\]{2Np) generators on mesons and baryons. Mesons consists of (P) and 1 (V) 

lowest-lying states, while baryons contain ^ (B) and | {B*) lowest- lying states. When this interaction is restricted 
to the sector PB PB it reproduces the standard WT off B targets. Its SF extension automatically yields the 
standard WT for PB* — >■ PB* (hence the compatibility between the two symmetries). Additionally, the extended 
interaction provides contact 5- wave vertices for VB VB, VB* VB*, PB o VB, PB o VB*, PB* o VB, 
PB* o VB* , and VB ^ VB*. As we have tried to argue above, these new vertices are well defined predictions 
of an approximate emergent symmetry of hadrons. So we adopt them as our starting point to describe interactions 
involving vector mesons. 

The Hamiltonian corresponding to the vertex in Eq. (2.8) can be written for any number of flavors and colors Nc 
[85]. For the physical case Nc = 3 

'^WT{^) = -Jp-{^^do^]^BBl,cDB''^'' A,B,... = 1,...,2Nf. (2.9) 

The indices A,B,..., denote spin and flavor, and so they take 2Nf values. ^^b{x) is the meson field, a 2Np x 2Np 
Hermitian matrix which contains the fields of 0~ (pseudoscalar) and 1~ (vector) mesons. This matrix is not traceless; 
for later convenience it includes the S\J{2Nf) singlet meson (the mathematical rji). The contribution of rji to $ is 
proportional to the identity matrix and so it does not couple in H^^t- The normalization of $(a;) is such that a mass 
term (with a common mass m for all mesons) would read im^tr($^). 

B{x) is the baryon field. B^^'^ is a completely symmetric tensor, that is, in the irrcp [3] of S\J{2Nf)- It has 56 
components for Np = 3, and 120 components for Np — 4, and contains the lowest-lying baryons with = i and 
I"*". The normalization of the field B is such that a mass term (with a common mass M for all baryons) would take 
the form Mj^B^^^B^^'^ . E.g. the fields B^'^^{x), B^^^{x)/V2, and B^^^{x)/V6 have the standard normalization of 
a fermionic field. We refer to the Appendix for the detailed construction of ^^b{x) and B^^'^{x) in terms of the 
individual meson and baryon fields for Np = 4. 

The Hamiltonian H^j. has precisely the same structure displayed in Eq. (2.6), this time for the SF group SU(2A^i?).'^ 

The predictions of the SF extended WT model for Np = 3 have been worked out in [87] for baryonic resonances, 
and in [110] for the mesonic version. Applications involving charm have been given in [81-83]. 

Before closing this subsection, we note that the eigenvalues of the relevant operator 4J|^ J|j in Eq. (2.8) can also 
be written using the quadratic Casimir operator of S\J(2Np) 

(4JmJb)m = 2(C2(M;2iVi.) - C2ifiM;^Np) - C2{pb;2Np)) . (2.10) 

For baryons in the irrep (with Young tableau) [3] of SU(2iVF) (56 or 120 for Np = 3 or 4, respectively) and mesons 
in [2,1^^^"^] (the adjoint representation of SU(2A''f), 35 or 63 for Np = 3 or 4, respectively), the baryon-meson 
states lie in the irreps [3], [2, 1], [5, P^^f-^] and [4, 2, 1'^^p-s-^^ which correspond to 56, 70, 700 and 1134 for Np = 3, 
and to 120, 168, 2520 and 4752 for Np ~ 4. The corresponding eigenvalues are [85] 

A[3] = -ANp, A[2,i] = -ANp - 6, A[5 i2iv_p-2] = 6, A[4_2,i2iVF-3] = -2. (2-11) 

The SF extended WT interaction is attractive in three multiplets and repulsive in the remaining one. The sum of all 
eigenvalues with their multiplicity, i.e. the trace of AJ\jJg, is zero, as follows e.g. from tr(J^j) = 0. The operator 
can be written as 

4JmJb = E^'^^'" (2.12) 

where ^ are the four baryon-meson irreps of SU(27Vi?), A^ the eigenvalues and the orthogonal projectors. This 
allows to compute the matrix elements using the SU(2iVi?) Clcbsch-Gordan coefficients [111]. 

C. Heavy-quark spin symmetry implementation 

Whereas the model just described can be used directly for three flavors, the extension to include charm requires 
more care. One reason is that chiral symmetry and SU(4) invariance are less reliable for fixing the interaction in the 



* Note that the singlet part in B\i-,j^B^^^ does not eouple since the matrix [$,c'o'J'] is traceless 
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sectors involving charm. At the same time, new specific symmetries of non rclativistic type arise when heavy quarks 
are involved [78-80]. In the heavy quark limit the number of charm quarks and the number of charm antiquarks are 
separately conserved. This implies a symmetry Uc(l) x Uc(l). Likewise, the terms in the QCD Hamiltonian which 
depend on the heavy quark or antiquark spin are suppressed, being of order l/rrih, where mh is the mass of the heavy 
quark. Therefore, in the heavy quark limit, arbitrary rotations of the spin carried by the c quarks and, independently, 
of the spin carried by the c antiquarks, would leave unchanged the energy of the hadronic state. ^ This implies a 
symmetry SUc(2) x SUc(2) in the heavy quark limit. These invariances are aspects of heavy-quark spin symmetry 
(HQSS). In what follows we refer to SUc(2) x SUe(2) x Uc(l) x Uc(l) as the HQSS group. 

The approximate HQSS reflects on the hadronic spectrum and for charm it has a level of accuracy similar to that of 
flavor SU(3). Taking HQSS into account implies that the model described in the previous subsection, the SF extended 
WT or just SU(8)-WT model (for four flavors), has to be slightly modified. In order to keep the model simple, we 
will impose exact HQSS on it. The alternative would be to introduce instead l/rrih suppressions in some amplitudes, 
but such improvement is beyond the scope of the present work. 

First, it should be noted that SF by itself already guarantees HQSS in many sectors. Consider for instance, the 
couplings involving the channels ND and ND*. These channels are related through HQSS since there should be 
invariancc under rotations of the c quark spin (leaving the light quarks unrotated), and this mixes D and D* . But 
the same invariance is already implied by SF, which requires symmetry under independent rotations of spin for each 
flavor separately. The only cases where SF does not by itself guarantee HQSS is when there are simultaneously c 
quarks and c antiquarks: SF implies invariance under equal rotations for c and c, but HQSS requires also invariance 
when the two spin rotations are different. 

To be more speciflc, let us consider baryon-meson channels, and let Nc be the number of c quarks and Nc 
the number of c antiquarks. Nc ranges from to 4, and iVg from to 1. SF guarantees HQSS in the sectors 
{Nc, Ns) = (0, 0), (0, 1), (1, 0), (2, 0), (3, 0), (4, 0), but not in the sectors (1, 1), (2, 1), (3, 1), (4, 1). As compared to the 
former sectors, the latter ones contain extra cc pairs. For the present discussion, we refer collectively to these sectors 
as sectors with "hidden charm", regardless of whether they have net charm or not. The hidden charm sectors are the 
main subject of the present work. We note that, for S'-wave interactions (the ones of interest here), even SU(6) SF, 
rather than SU(8), is sufficient to guarantee HQSS in the sectors without hidden charm: a rotation of the single heavy 
quark (or antiquark) can be produced by a light sector rotation followed by a global rotation, without changing the 
energy. In order words, in those sectors and for S'-wave. any SF invariant interaction enjoys HQSS automatically. 

It is perfectly possible to write down non trivial models enjoying simultaneously SU(8) and HQSS invariances 
(namely, by requiring SUg(8) x SUg(8)) but they would not reduce to WT in the light sector. Concretely SU(8)- 
WT conserves C = Nc ~ Nc but not Nc and TVg separately. Of course, one could impose this by hand, but it is 
automatically taken care of by our modifled interaction below (Eq. (2.22)). Also, the restrictions of SU(8)-WT to the 
sectors {Nc,Nc) = (1, 1), (2, 1), (3, 1), (4, 1) turn out to violate HQSS. 

In order to implement HQSS in the model let us analyze its content. We extract the trivial kinematic part and 
work directly in the space with only spin and flavor as degrees of freedom. Let 

ifwT =4Jl,J]j. (2.13) 

This operator can be written in terms of meson and baryon operators [81, 85], and it contains two distinct mechanisms 
which stem from expanding the meson commutator in Eq. (2.9), 

i^ac = - : A'/^^cM^bS^^^B]^^^ :, A, . . . , E ^ 1, . . . ,2Nf. (2.14) 

Here M^b and B^^'^ are the annihilation operators of mesons and baryons, respectively, with M^^b — {M^ a)\ 
and B^gfj = {B^^'~^y . B"^^^ is a completely symmetric tensor. They are normalized as 

[Af^B,A/tc^] ^ SA^c^ 

{B^^^,B\,g,^,} = 6i,S§,S^, + ■■■ (6 permutations). (2.15) 

Note that in the SU(8)-WT model, the rji (SU(8) singlet meson) does not couple and could be ignored, however, 
this meson has to be present in the corrected interaction since it mixes with the other mesons under HQSS. 



However all c quarks present in the state, being identical particles, are rotated by a common rotation, and similarly for the c. 
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FIG. 1: The two mechanisms acting in the spin-flavor extended WT interaction, il/ex (exchange of quarks) and H^c (annihilation 
and creation of quark-antiquark pairs). In the HQSS corrected version of the interaction, Eq. (2.21), the labels A and B in 7/ac 
only take light-flavor values. 



Schematically, representing the quark and antiquark operators by and Qa, 

M^B - Q^Qb. M^^b ~ Q^^Q^B^ B^""^ ^ Q^Q^'Q^, b\^c ~ QaQbQ^ ■ (2-16) 

So, an upper index in M or B represents the SF of a quark to be annihilated, whereas in A/^ it represents that of 
an antiquark to be created. Likewise, a lower index in or B'' represents the SF of a quark to be created while in 
M it represents that of an antiquark to be annihilated. From this identification it is immediate to interpret the two 
mechanisms i/ex and iJac in terms of quark and antiquark propagation. 

The two mechanisms involved, iJox and iJac sue displayed in Fig. 1. In Hex (exchange) the quark with spin- flavor 
A is transferred from the meson to the baryon, as is the quark with label B from the baryon to the meson. On the 
other hand, in iJac (annihilation-creation) an antiquark with spin-flavor B in the meson annihilates with a similar 
quark in the baryon, with subsequent creation of a quark and an antiquark with spin-flavor A. In both mechanisms 
the quarks or antiquarks C, D and E are spectators from the point of their spin-flavor (the ubiquitous gluons are 
not explicited). Also in both mechanisms effectively a meson is exchanged. In passing we note that the OZI rule 
is automatically fulfilled as regards to the exchanged meson. OZI rule violating mechanisms would be of the type 
depicted in Fig. 2 and are not present in WT. 




FIG. 2: OZI rule violating mechanisms. Gluons are implicit. 



It appears that 77ac can violate HQSS when the annihilation or creation of qq pairs involves heavy quarks, whereas 
HcK would not be in conflict with HQSS. This is indeed correct. To expose this fact more clearly, let us consider the 
symmetries of these two interaction mechanisms. Let Np ~ 4 and let [/ be a matrix of SF SU(8). Upper indices 
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transform with and lowers indices with U 

Qa W'aQb, Q\ ^ U^aQ^b- (2-17) 
Therefore (with obvious matrix/ tensor notation) 

M U^'MU, UHl^U, 

B ^ {U^ (g)U'' (E)U'')B, ^ B\U ®U ®U). (2.18) 

The indices are so contracted that iJox and _ffac arc both invariant under these SU(8) transformations. However, 
HQSS requires also invariance when the charm quark and the charm antiquark receive different rotations. To examine 
this, let us consider the transformation under U G SUq(8) and U £ SUq(8), i.e., different transformations for quarks 
and antiquarks (previously U = U). In this case 

Q^C/tg, Qt^i/tgt^ Q^QU, Q^^Q^U, (2.19) 

and therefore 

M UH'IU, A'/t -> UHI^U, 

B {U'' ®U'' (^U^)B, B^^B\U®U®U). (2.20) 

Clearly, the mechanism 77cx, which depends on the combination MM\ is still invariant under this larger group 
SUg(8) X SUq(8).^ It certainly preserves SF and HQSS. On the other hand, i/ac depends on the combination Af^M 
which transforms with f7, while BB"^ transforms with U . i/ac is SF invariant [U ~ U) but not HQSS invariant. A 
simple solution to enforce HQSS with minimal modifications is to remove just the offending contributions in i?ac, 
which come from creation or annihilation of charm quark-antiquark pairs. This implies to remove the interaction 
when the labels A oi B are of heavy type in i/ac-^ That is, we adopt the following modified version of the i?ac 
mechanism 

i/^, = - : A/t^cA/^^B^^^B^^^ :, C, D, E ^ I, . . . ,8, i,B = l,...,6. (2.21) 

The indices with hat are restricted to SU(6). By construction Nc and Nc are exactly conserved in H'^^. Also 
SUc(2) xSUc(2) is conserved: when U and U act only on the heavy sector M'^^cM^ ^ and ^^^^ are unchanged. 
So HQSS is preserved. On the other hand, when U = U and this matrix acts on the light sector, H'^^ is unchanged, so 
it enjoys SU(6) symmetry. Exact SF SU(8) and flavor SU(4) is no longer present. Presumably this breaking of SU(4) 
is comparable to the breaking through the kinematics due to the substantially heavier mass of the charmed quark. 
To summarize, our model (in all sectors) is given by 

K{s) I I I 

V = ^-^^WT> = Hex + H'^c- (2.22) 

This model fulfills some desirable requirements: i) It has symmetry SU(6) x HQSS, i.e., SF symmetry in the light 
sector and HQSS in the heavy sector, the two invariances being compatible, ii) It reduces to SU(6)-WT in the light 
sector, so it is consistent with chiral symmetry in that sector. And, iii) is a minimal modification that preserves 
simplicity and does not introduce new adjustable parameters. 

D. The model in the various sectors 

We can analyze the model in the different {N^ Nc) sectors, which as already noted, do not mix due to HQSS. 
In all the sectors without hidden charm, namely, {Nc, Ng) = (0, 1), (0, 0), (1, 0), (2, 0), (3, 0), (4, 0), H'^^ produces the 
same amplitudes as H^c when the latter is restricted to the corresponding sector. Indeed these interactions vanish 



® Note that the commutation relations Eq. (2.15) are also preserved by this symmetry. 
Keeping the contributions with A = B of heavy type would preserve Uc(l) x Uc(l), i.e., conservation of Nc and Nc, but not SUc(2) X 
SUc(2). 
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unless the state contains a quark-antiquark pair with quark and antiquark of the same type. In the absence of 
hidden charm, the pair must be hght and in this case the two operators produce the same result. This is consistent 
with our previous observation that, when there arc only heavy quarks or heavy antiquarks but not both, SF already 
implies HQSS. So in all these sectors, our model produces the same amplitudes as SU(8)-WT after removing channels 
involving hidden charm. This observation has been applied in [81-83, 88]. 

It is noteworthy that, in the sectors {Nc, Nc) = (0, 1) and (4, 0), corresponding to C = — 1 and C = 4, H!^^ = i/ac = 
as they do not contain light quark-antiquark pairs. Also these two sectors cannot couple to any other {Nc,Nc) sector 
in the baryon-meson case. Therefore for them, our model coincides directly with SU(8)-WT. 

Let us turn now to the sectors with hidden charm. These are {Nc,Nc) = (1, 1), (2, 1), (3, 1), (4, 1). For all these 
sectors H'^^ vanishes. The reason is that in these sectors the relevant quark-antiquark pair (quark and antiquark with 
equal labels) is necessarily of heavy type, and such amplitude has been removed from H'^^. (Note that i?ac does not 
vanish in these sectors.) So for the hidden charm sectors H^r^ reduces to the exchange mechanism H^x- 

The interaction is effectively i?ex for the four hidden charm sectors and also for C = — 1 and C = 4. This has the 
immediate consequence that the interaction i?-(vT ^^.s only two eigenvalues, namely —2 (attractive) and 6 (repulsive), 
in those sectors.^ This follows from the fact that Hex has a large (accidental) symmetry group, SUg(8) x SUg(8), 
which produces only two large multiplcts of degenerated states in the baryon-meson coupled-channels space. Under 
SUg(8) X SUg(8) the baryons fall in the irrep (120, 1) (120 being the symmetric representation of three quarks with 
8 possible spin-flavor states).^ Likewise, the mesons belong to (8,8*) (being qq states). Therefore, the baryon-meson 
states form two SUq(8) x SUq(8) multiplets: 

(120, 1) ® (8, 8*) = (330, 8*) e (630, 8*), (2.23) 

(330 is the symmetric representation [4] of SU(8) while 630 corresponds to the mixed symmetry [3,1]). The two 
corresponding eigenvalues are^° 

A(330,8*) = 6, A(630,8-') = —2. (2.24) 

These two eigenvalues are also present in the original SU(8)-WT interaction (namely, A2520 = 6 and A4752 = —2 from 
Eq. (2.11)) since the two interactions coincide for C = — 1 or C = 4. It can also be noted that under SF SU(8), these 
multiplets reduce as follows 

(330,8*) = 120® 2520, (630, 8*) = 120 ® 168 © 4752. (2.25) 

Of course, these are the same SU(8) irreps obtained from 120(8)63 (baryon-meson except 771) and 120 (g)l (baryon-?/i). 
Therefore, the interaction can be written using SU(8) Clebsch-Gordan coefficients by means of^^ 

Hex = 6(P{20 + ^2520) - 2(Pi20 + ^168 + ^4752)- (2.26) 

We emphasize that the large multiplcts (330,8*) and (630,8*) are not realized in our model. First of all, they 
contain sectors without hidden charm, for which the interaction does not reduce to the exchange mechanism H^x- 
And second, the eigenvalues 6 and —2 refer only to the driving operator AJljJg. The vertex V in Eq. (2.22) depends 
also on hadron masses and meson decay constants, for which we use physical values in V and in the propagators. The 
values that we use for these magnitudes are collected in Table II of [83] . 

E. Analysis of the hidden charm sectors 

Here we consider hidden charm sectors with C — 0,1, 2, 3, i.e. Nc ~ I and Nc ~ 1, 2, 3, 4, respectively. We want to 
classify the possible states under the symmetry group SU(6) x HQSS, with HQSS = SUc(2) x SUc(2) x Uc(l) x Uc(l). 
Since in the hidden charm sectors there is exactly one heavy antiquark, it is not necessary to specify the irrep of the 
factor SUc(2) x Uc(l) and we can use the notation R2,j^+i^c for the irreps of SU(6) x HQSS, R being the SU(6) irrep 
of the light sector, C the charm quantum number and Jc the total spin carried by the one or more c quarks (not 



Remarkably, these two eigenvalues are the only A'^^-independent ones in Eq. (2.11). 

Fermi-Dirac statistic is taken care of by the antisymmetric color wavefunction. 

They can be obtained by applying Hex to two suitable states, e.g., (Aft^iBjn i Mti2-B|j^j^)|0). 

One exception comes from the two 120 irreps, which differ by the type of symmetry of the four quarks. This information is not contained 
in the SU{8) Clebsch-Gordan coefficients. 
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including the spin of the c antiquarks). The corresponding dimension is dim il x (2Jc + 1) x 2 (the last factor coming 
from the two possible spin states of the c) . 

Subsequently, we study the breaking of light SF down to SU(3) x SU(2) keeping HQSS, and enumerate the number 
of attractive channels in each (C, r, J) sector, where r is the SU(3) irrep and J the total spin. 

In practice we will assume exact isospin and spin SU(2)/ x SU(2) j, as well as conservation of each flavor, but not 
exact SU(3) and HQSS, for the baryons and mesons composing the coupled-channels space. Therefore the sectors are 
labeled by (C, 5, /, J), S being the strangeness quantum number and / the isospin. This implies a further breaking 
of each (C, r, J) sector into (C, S, /, J) subsectors. 



1. C = 



For C = 0, the quark content is lUcc, with two possibilities of grouping into baryon- meson: {£M){cc) and {££c){£c). 
(Here £ denotes any light flavor quark, u, d, s.) The total dimension of the space is 56 x2x2 + 21x2x6x2 = 728, 
and contains the following SU(6) x HQSS multiplcts^^ 



n 



c=o 



56, 



I 56, 



The eigenvalues turn out to be 



A 



562 



A 



702 



70o 



= -2, 



(SU(6) X HQSS). 



A' 



= 6. 



(2.27) 



(2.28) 



The accidental degeneracy between 702, o and one 662,0 takes place in our model and it is not a necessary consequence 
of SU(6) X HQSS. This symmetry does not fix the three possible eigenvalues and the precise splitting between the 
two copies of 562,0- In general, the accidental degeneracy will be lifted in V and the T-matrix even when an exact 
SU(6) X HQSS invariance is assumed in the hadron masses and meson decay constants. 

Next, we consider the breaking of light SF SU(6) down to SU(3) x SUj,(2). E.g. 56 = 830104. HQSS is unbroken. 
After recoupling the spin carried by light and heavy quarks and antiquarks to yield the total spin J, we obtain 
representations of SU(3) x SU(2)j labeled as r2j+i, where r is the SU(3) irrep. This yields the following reductions 
(the two 562.0 have the same reduction) 



562,0 = (82 ' 
702,0 = (I2' 



' 104)2,0 = 
' 82 ® 84 



(82 ® 82 ' 
1 102)2,0 = 



584) 

(12 9 



3 (IO2 © 
I2 ® 14' 



1O4 © 1O4 g 

© (82 © 82 



lOe), 

B 84) © (82 



84©84©86)©(102©102l 



(2.29) 
IO4). 



In the reduction (82)2,0 = (82 © 82 © 84) in 662.0, the three octets only differ in how the light sector spin is coupled 
to the heavy sector spin, therefore these three irreps are degenerated if exact HQSS is assumed. (82 © 82 © 84) is a 
multiplet of SU(3) x HQSS. Similar statements hold in the other cases: each 662,0 produces two such multiplets and 
702,0 produces four. Consequently, in the hidden charm sector with C = we expect to find 8 different eigenvalues 
after SU(6) x HQSS is broken down to SU(3) x HQSS. Let Ai, A2 be the eigenvalues of the two multiplets in the 
repulsive 662,0, A3, A4 in the attractive 662.0, a-nd A5, Ag, Ay, As those in 702,o- In this case, the spectra in each 
(C, r, J) sector is as follows: 



12 


(A5,A5), 




14 


(As), 




82 


(Ai, Ai, A3, 


A3, Ag, Ag, Ay), 


84 


(Ai, A3, Ag, 


Ay, Ay), 


86 


(Ay), 




102 


(A2, A4, As, 


As), 


104 


(A2, A2, A4, 


A4, As), 


lOe 


(A2,A4). 





(2.30) 



In the SU(6) limit, Ai = A2, A3 = A4, A5 



Ay = Ag. Breaking down the symmetry to SU(3), one expects 



A 



3,4,5,6,7,8 



< < A1.2. 



(2.31) 



{M){cc) is purely 562,0 from the symmetry of the three hght quarks. The two hght quarks in {Uc) are symmetric giving a 21 of SU(6), 
coupled to the further light quark in (£5) gives 21 Cg) 6 = 56 © 70. These two 562,0 are not directly those in Eq. (2.28). 
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Each negative eigenvalue can give rise to a resonance or bound state. Each such state is a full multiplct of 
SU(3) X SU(2)j. This implies the following expected number of states in each (C, r, J) sector: up to two states in I2, 
one in I4, five in 82, four in 84, one in 8e, three in IO2, three in IO4, and one in lOe, all of them with C — 0. 

2. C = l 

For C = 1 there are two baryon- meson structures, namely, {£ic){cc) and {£cc){ic). The total dimension of the space 
is 384, with the following reduction under SU(6) x HQSS 



Hc=i - 21 



3,1 ' 



and eigenvalues 



2I31 ©2I1 



-^211,1 — AiE 



15 



3,1 



(SU(6) X HQSS), 
4i3. =6- 



(2.32) 
(2.33) 



Once again the accidental degeneracy beyond SU(6) x HQSS is lifted in the T-matrix. 
After the breaking SU(6) D SU(3) x SU(2)j^, and recoupling to J, one finds 



213,1 = (6393^)3,1 = (62062 
21i,i ^ (6303^)1,1 = (62064) 
153,1 = (61 3^)3,1 - (62064) 



64064066)0(3^034), 



(3^ 



(2.34) 



Thus for C = 1, each of the four SU(6) irreps split into two SU(3) x HQSS multiplets.^^ 

In principle there are eight different eigenvalues. Denoting by Ai,A2 the eigenvalues of the two multiplets in the 
repulsive 2X3,1, A3, A4 for the attractive 2X3,1, A5, A6 for 21i,i, and A7, Ag for I53.1, the spectra in each (C, r, J) sector 
is as follows: 



3; 


(A2, A4, 


Ae, As, As), 


3: 


(A2, A4, 


As, As), 


% 


(As), 




62 


(Ai, Ai, 


A3, A3, A5, A7), 


64 


(Ai, Ai, 


A3, A3, A5, A7), 


66 


(Ai,A3) 





(2.35) 

In the SU(6) limit, Ai = A2, A3 = A4, A5 = Ae, and A7 = As. If only SU(3) symmetry is assumed, we expect the 
following number of states: up to four states in 33, three in 3|, one in 3g, four in 62, four in 64, and one in 66. 



3. C = 2 



For the hidden charm sector with C = 2 there are two baryon-meson quark structures, {£cc){cc) and {ccc){£c). The 
space has dimension 120, with the following reduction and eigenvalues: 



-Hc- 



62.2 64,2 ( 
A64.2 = ~2, 



A' 



(SU(6) X HQSS) 
6. 



The SU(3) X HQSS muhiplcts are as follows: 



62,2 = (32)2,2 = (32 32 34), 

64,2 = (32)4,2 = (32 34 34 36). 



(2.36) 



(2.37) 



Note that the two (62 © 64) multiplets above do not mix if HQSS holds, as they carry different heavy quark spin labels. Such label has 
been left implicit to avoid clumsiness. 
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In this case the SU(6) x HQSS multiplcts are not reduced further under SU(3) x HQSS. Cahing Ai the eigenvalue 
of the repulsive 64 2, A2 that of the attractive 64.2, and A3 the one of 62.2, yields the following spectra for the three 
(C, r, J) sectors: 

82 : (Ai, A2, A3, A3), 
84 : (Ai, Ai, A2, A2, A3), 

36 : (Ai,A2). (2.38) 

This produces the following expected maximum number of states: three in 82, three in 84 and one in 85 

It can be noted in the present case, C = 2 with hidden charm, the assumption of SU(6) invariancc does not 
add anything (docs not reduce the number of parameters) on top of that of SU(3). The reason is that here HQSS 
automatically implies SF: there is just one light quark and rotations of it can be produced by combining global 
rotations with heavy quark rotations. 



C = 3 



For C = 3 there is just one quark structure, {ccc)cc. The dimension is 16. The SU(6) x HQSS reduction is 



Hc=3 - 13,3 ® 15,3 (SU(6) X HQSS), 



with eigenvalues 



The SU(3) X SU(2)j reduction is: 



A13 



13,3 = (11)3,3 = (I2© I4), 
15,3 = (11)5,3 = {UG)W- 

Once again, the SU(6) x HQSS multiplcts are not reduced further under SU(3) x HQSS. 
The following spectra is obtained for the various (C, r, J) sectors 



(2.39) 
(2.40) 

(2.41) 



12 


(A2), 


14 


(Ai,A2 


l6 


(Ai), 



(2.42) 

where Ai denotes the eigenvalue corresponding to 1.5,3, and A2 that of I3.3. So no states will be produced in Ig, and 
up to one state is expected in I2 and I4. 

Some of the results of this subsection arc summarized in Table I. 





C 





1 


2 


3 


J 


SU(3) 


1 8 10 


3* 6 


3 


1 


1/2 




2 7 4 


5 6 


4 


1 




(2) (5) (3) 


(4) (4) 


(3) 


(1) 


3/2 




1 5 5 


4 6 


5 


2 




(1) (4) (3) 


(3) (4) 


(3) 


(1) 


5/2 




1 2 
(0) (1) (1) 


1 2 
(1) (1) 


2 

(1) 


1 

(0) 



TABLE I: Total number of channels for each J and each SU(3) irrep, for the various hidden charm sectors. Here each channel 
represents a full SU(3) and spin multiplet (C, r, J) (rather than a isospin-spin multiplet, {C, S, I , J).) The expected number 
of resonances in each sector is shown between parenthesis. The actual number of resonances will depend on the values of the 
physical masses and meson decay constants. 
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F. Lagrangian form of the interaction 



111 the strict heavy quark Hmit, the total spin J, and the separate spins of the hght, heavy quarks and heavy 
antiquarks subsystems are conserved, when only S'-wave interactions are considered. Moreover, the matrix elements 
of QCD Hamiltonian depend only on the spin-flavor quantum numbers of the light degrees of freedom. Each n-fold 
degenerated SU(6) or SU(3) multiplet implies the specification of n{n + l)/2 parameters which are coupling constants 
of the corresponding operators present in the interaction. From the previous analysis it follows that the number of 
independent operators in the hidden charm sectors with C = 0,1,2,3 is 4, 5, 4, 2, respectively, for a generic interaction 
if SU(6) X HQSS invariance is assumed and 12, 10, 4, 2 if the symmetry is reduced to SU(3) x HQSS. 

In what follows we will focus on the sector with hidden charm and C = 0. For this sector we will write down the most 
general (modulo kinematical factors) S'-wave Lagrangian consistent with SU(3) x HQSS for the baryon-meson coupled- 
channels space. This Lagrangian contains 12 operators and our model gives well defined values for the corresponding 
coupling constants. We should note that we actually compute the matrix elements of our interaction using directly 
the previous expressions, either in terms of projectors using Clcbsch-Gordan coefficients, or of hadron creation and 
annihilation operators in spin-flavor space, by taking Wick contractions. Nevertheless, writing the interaction in 
field-theoretical Lagrangian form is highly interesting in order to make contact with alternative approaches in the 
literature. 

For this purpose it is convenient to organize the hadroiis forming multiplets of HQSS into building blocks with well 
defined HQSS transformation properties [80, 112-114]. Specifically, consider a HQSS doublet composed of pseudoscalar 

meson and vector meson with one heavy quark (e.g. D and D*). Let M'^'^'> and Mjf'' be the corresponding fields, then 

M(^) = 7oM('=)t7o = (A/(")J,(-)7^ - A'/(")t(-)^5)g_^. (2.43) 

As usual (±) represent the positive and negative frequency part of the fields, corresponding to purely annihilation for 
(-I-) and purely creation for (— ). Therefore, Ai'^^^ [Af^^^] annihilates [creates] the meson with one heavy quark but it 
does not create [annihilate] the corresponding antimeson with a heavy antiquark. A similar proviso is applied in all 
the other fields for hadrons carrying heavy quarks and/or antiquarks. Besides, 

1± 1^ 

Q± = ^ (2.44) 

where v'^ is the heavy hadron velocity (u^ = 1). We use Bjorken and Drell [115] conventions for the Dirac gammas. 
For the hadron fields we use the conventions of [111] and this fixes the relative sign between pseudoscalar and vector 
(see Appendix A). 

Likewise, for a HQSS meson doublet with one heavy antiquark 

M(^) = (A/^"")(+)7^-HA/^"'(+)75)Q_, 

= 7oM(^)^7o = Q-(M('^'t (-'7'^ -M('^)^(-)75). (2.45) 

For a HQSS meson doublet with one heavy quark and one heavy antiquark (e.g., r]c and J/ip)'- 

M("^) = g+(A4'=^)(+)7'^ + A/('="")(+)75)Q_, 

M("^) = 7oM("^)t7o =Q_(Af('='^)|/-)7^- Af(==)'f(-)75)g+. (2.46) 
For a HQSS baryon doublet with exactly one heavy quark (e.g., Ec and E*), 

v3 

v3 

Here B^'^' is the Dirac spinor of the 1/2+ baryon in the doublet while _b('=)^ is the Rarita-Schwinger field for the 3/2+ 
baryon: u^B^'^)^ = ^^,3^"'^'' = 0. 

Finally, for a HQSS singlet baryon with exactly one heavy quark (e.g., Ac) 

= B'-"^ = B(")t^o = ^^"^^"^ (2.48) 
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In addition, the polarization of the baryons carrying heavy quarks is such that (5_i?('^)(+) = (5_i3('^^^(+) = 0, hence 

Q_B(^) = 0_B(='^ = 0. (2.49) 
Under HQSS rotations these fields transform as follows 

ScM'^'^K M^'^^S-j-i, ScM^"""^ S^^ , ScB^^l ScB^"'"' 

M^'''>S-\ ScM^^'l S-Ja'^"''^ S-^ , B("'S'-\ B^'^^S--!. (2.50) 

Here Sc and Sc are the matrices in Dirac space representing the c or c spin rotation and satisfy S'^^g = 70'S'c.c7o- 
All other hadrons without heavy quarks nor antiquarks are HQSS singlets. They have complete fields (positive and 

negative frequency parts) and are denoted without boldface type. 

The hadrons are also organized into SU(3) multiplets. We use fields with labels a,b,c,... = 1,2,3 (or up, down, 

strange) in the fundamental or anti- fundamental representations of SU(3), in such a way that 

T^--- ^u^^a-W^'b-'-T^,':::, uesu{3) (2.51) 

For C ~ with hidden charm, the following SU(3) multiplets are needed 

£,a ^ (^^0 (2.52) 

Here represents the HQSS doublet formed by and _D*°, etc. 

S^"^ = V2£^° (2.53) 

This is a symmetric tensor (irrep 6 of SU(3)). S^, H^, and $7^ arc the HQSS doublets with (Sc,S*), (S^,S*), and 
(ricjf^c), respectively. 

= (-H^ H+ -A,) (2.54) 

In this case the members of the SU(3) multiplct are HQSS singlets. Further we define as the SU(3) singlet and 
HQSS doublet containing {r]c,J/ip). 

In addition, the following light baryon multiplets appear: with = 1/2+ 



V6 V2 

■=■0 /2 , 



(2.55) 



and with = 3/2+, A"*"^, a Rarita-Schwinger field and a symmetric tensor normalized as 



Al" = ^/6A++, Aii2 = V2A+, A;,^^ = V2A0, Af ^ = V6A;, 
All3 = ^/2S;+, Ai23^i];o^ A223^y2S;-, 

Ai33^V2s;°, A^33 = y2s;-, 

A333 = (2.56) 

The relative phases of all fields here arc standard with respect to the conventions adopted in [111] for the rotation, 
flavor and spin-flavor groups. So for instance, (S+,S]?) is a standard isospin doublet and (I]+,S'',S^) is a standard 
isovector. For SU(3) (and SU(4)) the convention in [116] is used instead of that in [117].^^ 



^'^ The matrix elements between standard states of the step operators « <-> d, d s, and s •<-> c arc required to be non negative [116], 
rather than those of u <r> d and u o s [117]. 
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Regarding parity, wc note that the hadrons with spin-parity 1 or 1/2"*" have normal parity, whereas those with 

0~ or 3/2+ have abnormal parity. So S''^, V'l ^ca, and 5° are true tensors while A"'"^, and S^afa are pseudotensors. 

Also, T,^a, Sea, A"**^, and S^ab have large upper components while and 5° have large off-diagonal blocks in Dirac 

space. The 75 matrices introduced to preserve parity fit in this scheme. 

Next, we write down the 12 most general operators allowed by SU(3) x HQSS in the baryon-meson coupled-channels 

•«-> — > <~ 

space, in S'-wave and preserving parity. The operator dy = w^(i9^ — 9^) acts on the mesons only and it is introduced 
in order to produce the correct kinematical dependence in the amplitudes. 



C2{x 

£3(2; 
Li(x 
Cz{x 

Ce{x 
L-j{x 

Cg{x) + Cio{x 
Cii{x) +Ci2ix 



h.c, 



.giS^E^tr(^za.i/»), 

g3Hc^Ha„)5bI]\+h.c., 

^ -h.c, 

gs e'-^^i^abS^d tr(5e7^75 iduD'^) + h.c. 



55 is^,b-(/'HE)£)c A; 



1, 



^'iabK tr(Dc(g9 <?M^ + .910 ia^,^)idyD''), 



S^acSf" tr(£>b(gn 5m- + 512 ia^,) id^D"). 



The traces refer to Dirac space. 

The reduction of these Lagrangians when no strangeness is involved is as follows: 



Ci{x) 




C2{x) 


= 52 A^A^tr(^iU,i/'), 




= g3A^xP{-id.,)DN + h.c., 


£4(2;) 


= g4'E'^jip{-idy)DTjjf,j5N + h.c., 


£5(2;) 


- S^i^(-i9,)-0r_,-A^ + h.c. 


£6(2;) 


= geAcActiiDidyD), 


£7(2^) 


= 0, 




<-> 

= .98 Ac tT{D Tjj^j^ i d„D) + h.c 



£9(2;) + £10(2;) + £11(2;) + Ci2{x 



(2.57) 
(2.58) 

(2.59) 
(2.60) 
(2.61) 

(2.62) 
(2.63) 
(2.64) 
(2.65) 

(2.66) 



(2.67) 

(2.68) 

(2.69) 

(2.70) 

(2.71) 

(2.72) 
(2.73) 

(2.74) 



S^jS^, tr(r>(G9 .9p. + Gio ia^,)idyD) 
-HS^jS^fc tr{DTjTk{Gii 9t^'^ + G12 i(T^,) idyD) . (2.75) 
Here j, k are isovector indices, t are the Pauh matrices and (3/2, A/lexf t|l/2, m) = C(l/2, 1, 3/2; m. A, M). Further, 



G9 — .99 4- 2911, Gio — gio + 2512, Gil — G12 — —gi2- 

Our WT model with SF and HQSS gives the following values for the parameters: 



(2.76) 



91 = 0, .92 = 0, .93 



2' 



54 



57 = -^, .98 = ^, .99 = 0, 510 = 0, 



55 

.911 = 



-1, 56 



1 



1 _ 1 

2' ^''-^2- 



(2.77) 



where we have defined 5^ = ^f^gi- 

The vanishing of gi and 52 follows from the OZI rule, which is fulfilled by the model. 
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III. COUPLED-CHANNELS UNITARIZATION AND SYMMETRY BREAKING 



A. Unitarization and renormalization scheme 



As previously discussed, the baryon-meson interaction is mediated by the extended WT interaction of Eq. (2.22) 
that fulfills SU(6) x HQSS and it is consistent with chiral symmetry in the light sector. The final expression for the 
potential to be used throughout this work is 

^r" = Af"4^(fc."+fc;°), (3.1) 

where and fc^" arc the CM energies of the incoming and outgoing mesons, respectively, and fi and fj are the 
decay constants of the meson in the i-channel and j-channel.^^We use the hadron masses and meson decay constants 
compiled in Table II of Ref. [83] . In particular, /j/* is taken from the width of the J/^* e~e^ decay, that is, fj/^ — 
290 MeV and we set /,,^ = /,//*, as predicted by HQSS and corroborated in the lattice evaluation of Ref. [118]. The 
Dij are the matrix elements of i?(vTj (2-22), for the various hidden charm CSIJ sectors previously discussed. 
Those for the strangeless hidden charm C = case, for which i^^^T = Hex, and that we will discuss in what follows 
are given in Appendix B. 

In order to calculate the scattering amplitudes, Tij, we solve the on-shcll Bethc-Salpcter equation (BSE), using the 
matrix v'~^^^'^ as kernel: 

rpCSIJ — (1 _ yCSIJ QCSIJ-^-lyCSIJ ^3 2) 

where G'-^^^"' is a diagonal matrix containing the baryon-meson propagator for each channel. Explicitly 

Gg^"(,s) = (Jo(V^;Af„m,) - Jo(/i^^; Af„ m,)) , (3.3) 

Mj {rrii) is the mass of the baryon (meson) in the channel i. The loop function Jq can be found in the appendix of 
[48] (Eq. A9) for the different possible Riemann sheets. The baryon-meson propagator is logarithmically ultraviolet 
divergent, thus, the loop needed to be renormalized. This has been done by a subtraction point regularization such 
that 

Gg^^-^(,s) = at V~s^^^'"'', (3.4) 



with /i*-^"^^ ~ \J m^j^ + M'^Yi 1 where TOth and Mth, are, respectively, the masses of the meson and baryon producing the 
lowest threshold (minimal value of mth + -^^^th) for each CSI sector, independent of the angular momentum J. This 
renormalization scheme was first proposed in Refs. [65, 66] and it was successfully used in Rcfs. [56, 81, 83]. A recent 
discussion on the regularization method can be found in Ref. [60]. 

The dynamically-generated baryon resonances appear as poles of the scattering amplitudes on the complex energy 
•v/i plane. One has to check both first and second Riemann sheets. The poles of the scattering amplitude on the first 
Riemann sheet that appear on the real axis below threshold are interpreted as bound states. The poles that are found 
on the second Riemann sheet below the real axis and above threshold are identified with resonances^^. The mass and 



As compared to our previous work of Ref. [83], we have 

• approximated (2-^/5 — Mi — Mj), with Mi and Mj the incoming and outgoing baryon masses, by the sum of the CM energies of 
the incoming and outgoing mesons. In the present case, the non relativistic approximation tends to increase the binding by up 
to few tens of MeV. This non-relativistic approximation for the baryons is consistent with the treatment for the baryons adopted 
in the previous section to implement the HQSS constraints, and it makes easier to connect with the effective HQSS Lagrangians 
introduced in Eqs. (2.57)-(2.66). 

• and, also for this latter reason, moved the + M)/{2M) factors included in the potential used in [83] to the definition of the 
loop function in Eq. (3.3). 

For the sake of completeness, and to make possible the determination (sec Eq. (2.77)) of the coupling 37 of the HQSS effective Lagrangian 
of Eq. (2.63), we also give in the Appendix B the coefficients for the rest of hidden charm sectors with explicit strangeness. 
Often we refer to all poles generically as resonances, regardless of their concrete nature, since usually they can decay through other 
channels not included in the model space. 
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the width of the state can be found from the position of the pole on the complex energy plane. Close to the pole, the 
scattering amplitude behaves as 

T-^-\s)^^-^p^. (3.5) 

The mass AIn and width Ffl of the resonance result from ^/sB, = A/j^ — irp{/2, while gje^"^^ (modulus and phase) is 
the coupling of the resonance to the j-channel. 



B. Symmetry breaking 

As it was already pointed out in the subsection HE, we classify states under the symmetry group SU(6) x HQSS, 
and consider the breaking of the light SF SU(6) to SU(3) x SUj, (2). Subsequently we break the SU(3) light flavor 
group to SU(2) isospin symmetry group, preserving the HQSS, and finally we break the HQSS. Thus, we assume 
exact isospin, total spin and flavor conservation. The symmetry breaking is performed by adiabatic change of hadron 
masses and meson weak decay constants, as it was previously done in Ref. [83]. At each symmetric point, the hadron 
masses and meson decay constants are averaged over the corresponding group multiplets. Further, we introduce three 
parameters, x,x' and x" that are changed from to 1, to gradually break the symmetry from SU(6) x HQSS down 
to SU(3) X HQSS, then to SU(2) x HQSS, and finally down to SU(2) isospin, respectively: 

m{x) = (1 - x) mgu(6)xHQSS + X msu(3)xHQSS, 
f{x) = (1 — x) /sU(6)xHQSS + X /sU(3)xHQSSj 
m(x') = (1 - x') TOsU(3)xHQSS + »71sU(2)xHQSS, 
f{x') = (1 - x') /sU(3)xHQSS + x' /sU(2)xHQSS, 
■m{x") = (1 - x") TOsU(2)xHQSS + x" TOSU(2), 

fix") ^{1- x") /su(2)xHQSS + X" /su(2). (3.6) 

In this way we can assign SU(3) and SU(6) representation labels to each found resonance, and also identify the HQSS 
multiplets. We will show below a diagram (Fig. 3) with the evolution of the hidden charm N and A pole positions as 
the various symmetries are gradually broken. 



IV. CHARMLESS AND STRANGELESS HIDDEN CHARM SECTOR: THE N AND A STATES 

In this work we will only discuss results on hidden charm baryon resonances with total charm C = and strangeness 
S ~ Q. Other sectors with charm different from zero will be studied elsewhere. 

In this sector, we find several 7=1/2 and / = 3/2 states, which correspond to A^-like and A-like states, respectively 
(here we use the same notation as in Refs. [74, 75]). All these states have odd parity and different values (J = 1/2, 3/2 
and 5/2) of total angular momentum. The list of coupled-channels and the corresponding coefficients Df^ can be 
found in the first six tables of Appendix B. 

In this hidden charm sector and in the SU(6) x HQSS limit, we saw (Eqs. (2.27) and (2.28)) that the group structure 
of the HQSS-constrained extension of the WT interaction developed in this work consists of two 562,0 and one 702,o 
representations. One of the 562,0 multiplets and the 702,o one are attractive. Thus, from the decomposition in 
Eq. (2.29) (see also Table I), we could expect up to a total of ten iV-like and seven A-like resonances.^® Because of 
the breaking of the SU(6) x HQSS symmetry, due to the use of physical hadron masses and meson decay constants, 
we only find seven heavy N and five heavy A states in the physical Riemann sheets. They have masses around 
4 GeV and most of them turn out to be bound. The remaining missing states show up in unphysical Riemann sheets. 
The evolution of all states as we gradually break the symmetry from SU(6) x HQSS down to SU(3) x HQSS, then 
to SU(2) X HQSS, and finally down to SU(2) isospin, is depicted in Fig. 3. Thanks to this latter study, we could 
assign SU(6) x HQSS and SU(3) x HQSS labels to each of the predicted resonances, which are all of them collected 
in Tables II and IV, and could also identify two HQSS multiplets in each isospin sector. 



Those lie in the SU(3) octets and decuplets irreps, respectively, contained in the attractive 662.0 and 702,o multiplets. 
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A. AT states (C = 0, 5 = 0, I = 1/2 ) 

As mentioned above, the model predicts the existence of seven heavy nucleon resonances: three states with the 
spin-parity = ^ , also three states with | sectors, and one state with = | . Their masses, widths and 
couplings to the different channels are compiled in Table II. 

• J ~ 1/2: In this sector, there are seven coupled channels, with the following threshold energies (in MeV): 

Nr]c NJ/^j AcD AcD* T.cD ^cD* 
3918.6 4035.8 4153.7 4294.8 4320.8 4461.9 4526.3 

• J = 3/2: In this sector, there are five coupled channels, with the following threshold energies: 

NJ/ij AcD* S*Z) J^cD* S*Z)* 
4035.8 4294.8 4385.2 4461.9 4526.3 

• J = 5/2: In this sector there is only one channel, I]*Z)*. with threshold equal to 4526. 3MeV. 

From the group decomposition of the SU(6) x HQSS representations, we could expect up to a maximum of five 
states with spin J = 1/2 (see Table I): one state from each of the two J = 1/2 octets encoded in the attractive 
562,0 representation, and three states corresponding to the each of the 82 octets that appear in the reduction of the 
702,0 representation [Eq. (2.29)]. However, the two poles related to the 562,0 representation appear in an unphysical 
Ricmann sheet, at the physical point (i.e., at the point of the evolution when the hadron masses and meson decay 
constants attain their physical values). As it can be seen from Fig. 3, these poles dissappear from the physical 
sheet when we pass from the SU(3) x HQSS limit to the SU(2) x HQSS one. Indeed, we could observe how the 
(82)2,0 C (562)2,0 pole almost coincides with the threshold value of the degenerated Nrjc and NJ /ip channels in the 
first steps of this evolution until it finally disappears. On the other hand, the (82)2,0 C (562)2,0 pole also gives rise 
to an octet of J = 3/2 states (see Eq. (2.29)), which is also lost at the physical point. Thus, for J = 3/2 we are also 
left only with the three baryon resonances stemming from the 702,o representation, one from (82)2.0, and two from 
(84)2.0- The J = 5/2 state is originated also from this latter multiplet. 

From the above discussion, it is clear that the A^-like resonances found in this work, and collected in Table II, form 
two HQSS multiplets. In the first one the light degrees of freedom have quantum numbers (82)2.0 C (702)2,0- This 
multiplet is formed by the three first resonances of the table (two with spin 1/2 and third one with spin 3/2) that 
correspond to the blue states in Fig. 3. They only differ in how the light sector spin is coupled to the spin of the 
cc pair. The second HQSS multiplet corresponds to (84)2,0 C (702)2,0 quantum numbers for the light sector, and it 
consists of the four remaining states in Table II (displayed in green in Fig. 3) one with spin 1/2, two with spin 3/2, 
and another one with spin 5/2. 

The members of each HQSS multiplet are nearly degenerate, but not totally because we also break the HQSS by 
the use of physical hadron masses. 

A word of caution is needed here. The mass of the J — h/2 resonance is around 4027.2 MeV. In this sector there 
is only one channel (E*i5*, with threshold equal to 4526.3 MeV), thus this state is around five hundred MeV bound. 
We expect our model to work well close to threshold, and therefore, in this case, interaction mechanisms neglected 
here and involving higher partial waves could be relevant for determining the actual properties of this resonance. 

There exist previous works on hidden-charm odd-parity nucleon states, also named crypto-exotic hadronic states. 
These studies can be divided in two types. Namely, those based on a constituent quark description of the resonances, 
and those where they are described as baryon-meson bound molecules or resonating states. Some of the predictions 
of these other models are compiled in Table HI. 

The baryon-meson coupled-channel calculations by Hofmann and Lutz for = 1/2^ in Ref. [65] and for = 3/2~ 
in Ref. [66] rely on a model of zero-range t-channel exchange of light vector mesons, based on chiral and large Nc 
considerations, and supplemented with SU(4) input in some vertices. This model is used as the driving interaction of 
pseudo-scalar mesons with the = 1/2^, 3/2+ baryon ground states. After solving the BSE using a renormalization 
scheme similar to that proposed here, some l/2~, 3/ 2 "-resonances were dynamically generated. Vector mesons in the 
coupled-channel space were omitted in those early studies, thus channels like TicD* or AcD* were not considered. 

More recently baryon-meson calculations using a hidden-gauge model have been carried out in Refs. [74-76]. These 
works consider 1 /2+ baryons interacting with pseudoscalar mesons and dynamically generate — hidden-charm 
nucleon resonances as poles in the T-matrix. Yet, the interaction of vector mesons with 1/2+ baryons {VB — >■ VB^ 
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SU(6) X HQSS 
irrep 


SU(3) X HQSS 
irrcp 


Mfl[MeV] 


rfl[MeV] 


Couplings 
to main channels 


J 


702,0 


(82)2,0 


3918.3 


0.0 


gNva = 0.5, gNj/i, = 0.6, gjy^o = 3.1, 3a,S' = 0.5, 
5e,S = 0.2 , g^^ = 2.6, g^* s* = 2.6 


1/2 


702,0 


(82)2,0 


3926.0 


0.1 


gN,,e = 0.2, gNj/i, = 0.04, Pa,S = 0.4, 5a,S* = 3.0, 
te,D = 4.2, = 0.2, = 0.7 


1/2 


702,0 


(82)2,0 


3946.1 


0. 


gNj/i, = 0.2, 3a,d' = 3.4, flsjD = 3.6, Se.d* = 1-1, 
= 1-5 


3/2 


702,0 


(84)2,0 


3974.3 


2.8 


gNrje = 0.5, 5jvj/v< ~ 0.05, gAe-S = 0.4, Sa.S* = 2.2, 
5e,d = 2.1, 3e,S' = 3.4, gs;B* = 3.1 


1/2 


702,0 


(84)2,0 


3986.5 


0. 


gjv.//v< = 0.2, 3a,S' = 1-0, 5ejS = 2.7, 5e<,S* = 4.3, 
= 1-8 


3/2 


702,0 


(84)2,0 


4005.8 


0. 


5JVJ/V. = 0.3, 3a,d' = 1-0, 5ejo = 1-6, 5e<,d* = 3.2, 
g^.B, = 4.2 


3/2 


702,0 


(84)2,0 


4027.1 


0. 


Pe; d' = 5.6 


5/2 



TABLE II: Odd parity hidden charm N {J = 1/2, J = 3/2 and J = 5/2) resonances found in this work. The first two columns 
contain the SU(6) x HQSS and SU(3) x HQSS quantum numbers of each state, while Mr and Tr stand for its mass and width 
(in MeV). The largest couplings of each pole, ordered by their threshold energies, are collected in the next column. In boldface, 
we highlight the channels which are open for decay. Finally, the spin of the state is given in the last column. Resonances with 
equal SU(6) x HQSS and SU(3) x HQSS labels form HQSS multiplets, and they are collected in consecutive rows. 

is also taken into account in [74-76], which leads to additional and degenerate — 1/2^ and 3/2^ hidden-charm 
nucleons. However, the J = 3/2+ baryons are not included at all, and thus some channels like are excluded. 

The main difference among our scheme and the hidden-gauge models is the definition of the coupled-channels 
space. We consider simultaneously pseudoscalar meson-baryon {PB) and vector meson-baryon {V B) channels, with 
= 1/2+ and 3/2+ baryons. However in the approaches of Refs. [74-76] all interaction terms of the type PB — >■ VB 
are neglected. Furthermore, channels with = 3/2+ baryons are not considered cither. The potential used in 
[74-76] for the PB — > PB transitions, with = 1/2+ baryons, is similar to that derived here. However, there 
exist important differences in all transitions involving vector mesons. When restricting our model to the PB — > PB 
sector, we still do not obtain the same results as in Refs. [74-76]. This is mainly due to i) the use of a different 
renormalization scheme and, ii) the presence in these latter works of a suppression factor in those transitions that 
involve a i-channel exchange of a heavy charm vector meson. 

However, when we use our full space, the inclusion of a similar suppression factor in our HQSS kernel is not 
quantitatively relevant for the dynamical generation of the resonances. Note that HQSS does not require the presence 
of such suppression factor. In summary, when we compare our approach with the other molecular-type ones, we 
observe in our model a rich structure of resonances due to the many channels cooperating to create them. Our states 
are much more lighter that those predicted in Refs. [74-76], though significantly less bound that the crypto-exotic 
baryons reported in Refs. [65, 66]. 

Finally, we will pay attention to the recent work of Ref. [119]. There a constituent quark model is used to describe 
isospin 7 = 1/2 baryons with uudcc quark content. The mass spectra is evaluated with three types of hyperfine 
interactions: color-magnetic interaction (CM) based on onc-gluon exchange, chiral interaction (FS) based on meson 
exchange, and instanton-induccd interaction (INST) based on the non-perturbativc QCD vacuum structure. The FS 
(CM) model predicts the lowest (highest) mass for each state. Results for the FS and CM models are displayed in 
Table HI. In all cases, the mass predicted by the INST model (not displayed in the table) lies between the values 
predicted by the other two models. Our results are closer to those predicted by the FS model, specially for the lowest 
lying states. 

B. A states (C = 0, 5 = 0, 7 = 3/2 ) 

The model predicts in this sector the existence of five heavy resonances (bound states; all of them appear below 
threshold): three with spin-parity = ^ and another two with = | . Their masses, widths and couplings to 
the different channels are compiled in Table IV. 
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M |MeV| 



\ KD_, 
\ 



4200 - 



4000 



SU (6) xHQSS ^ 
SU(3)xHQSS 



Sl^(3) xHQSS -> 
S^(2)xHQSS 







(102)2,0 


\ NJ/>// 




L \ 


(84)2,0 


V \ 




\ \ 






(82)2,0 


\ w* 




\ w 




SU(2)xHQS^^ 



SU (2) 



A_,,2 (4306.8) 
Al/2 (4306.2) 



Al,2 (4050.0) 
A3,2 (4032.5) 
)V5/2 (4027.2) 
Ai,, (4005.8) 
Nii2 (4005.8) 
N,i2 (3986.6) 
W,/2 (3974.3) 

^3/2 (3946.1) 
^1/2 (3926.0) 
Win (3918.3) 



^> (104)2,oc562,o, J=3/2, weakly bound 

^» (104)2,oc562,o, J=l/2, weakly bound 

A, (102)2,oc702,o,J=3/2 

A, (102)2,oc702,o,J=l/2 

N,(84)2,oc702,o,J=5/2 

N, (84)2,oc702,o,J=3/2 

N,(82)2,oc702,o,J=3/2 

N,(84)2,oc702,o,J=l/2 

N,(82)2,oc702,o,J=l/2 

Thresholds 

(104)2,oc562,o , weakly bound 

(102)2,0C702,0 

(82)2,oc562,o, weakly bound 

(84)2,0C702,0 
(82)2,0C7020 



FIG. 3: Evolution of the poles as symmetries, starting from SU(6) x HQSS, are sequentially broken to reach the isospin 
symmetric final crypto-exotic A'^ and A odd-parity resonances. The meanings of x, x' and x" can be found in (3.6). The lower 
index of the final states stands for the spin J of the corresponding resonance. The thresholds (red dashed lines) are marked 
together with the respective baryon-meson channel. The SU(6) x HQSS labels 702,o and 662,0, and the SU(3) x HQSS labels 
(82)2,0, (84)2,0, (102)2,0, (104)2,0 are also shown at the corresponding symmetric points. 



• J — 1/2: In this sector, there are four coupled channels, with the following threshold energies (in MeV): 

AJ/^/> EcD Eci5* E*i5* 
4306.9 4320.8 4461.9 4526.3 



• J = 3/2: In this sector, there are five coupled channels, with the following threshold energies: 



A77c AJ/V' ^ID EcL** E*^* 
4189.7 4306.9 4385.2 4461.9 4526.3 

• J = 5/2: In this sector there arc only two channels, with the following threshold energies: 



AJ/i/; E*L)* 
4306.9 4526.3 
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Ref. Model 


N{l/2~) 

[MeV] KcD AcD* EcD EcD* E*i3* 


iV(3/2-) 

1 v± /-^ ly vjvj ±i±cl±i± Vjiicumido 

[MeV] AcD* E*i3 EcD* E*7?* 


iV(5/2-) 
[MeV] T,*D* 


This work (82)2,0 C (702)2,0 
(84)2,0 c (702)2,0 


3918 3.1 0.5 0.2 2.6 2.6 
3926 0.4 3.0 4.2 0.2 0.7 


3946 3.4 3.6 1.1 1.5 




3974 0.4 2.2 2.1 3.4 3.1 


3987 1.0 2.7 4.3 1.8 
4006 1.0 1.6 3.2 4.2 


4027 5.6 


[65] [66] zero-range 

vector exchange 


3520 5.3 


3430 5.6 




[75] hidden-gauge 


4265 0.1 3.0 
4415 0.1 2.8 


4415 0.1 2.8 




[76] hidden-gauge 


4315 X X 
4454 X X 


4454 X X 




[119] quark model 
uudcc 


FS-CM 

3933-4267 

4013-4363 

4119-4377 

4136-4471 

4156-4541 


FS-CM 

4013-4389 
4119-4445 
4136-4476 
4236-4526 


4236-4616 



TABLE III: Comparison of the whole spectrum of hidden-charm nucleons (or crypto-exotic nucleons) with odd-parity and 
angular momentum L = predicted by our model with some results from previous models. In all cases, masses and couplings 
(g) to the dominant channels (when available) are shown in sequential rows. In the hidden gauge model of Ref. [76], the 
numerical values of the couplings are not given. In this case, we indicate with a symbol X the elements of the coupled channel 
space used to generate each resonance. On the other hand, in the case of the predictions of this work, in the column "Model" 
we give the HQSS multiplet. Besides, we have also omitted the small couplings to the Nrjc and NJ/ip channels that can be 
seen in Table II. 

We obtain three A(J = 1/2) states as expected from the group decomposition of the SU(6) x HQSS representations 
(see Table I): one state from each of the two J = 1/2 decuplets encoded in the attractive 702,o representation, 
and a further state corresponding to the J ~ 1/2 decuplet that appears in the reduction of the 562.0 representation 
[Eq. (2.29)]. The evolution of the corresponding poles is shown in Fig. 3. 

The pole that corresponds to (104)2,0 C 662,0 (light magenta circles in Fig. 3) has a mass quite close to the Arjc 
and AJ/ip degenerated thresholds, between the SU(6) x HQSS and the SU(2) x HQSS symmetric points. Later, 
while moving to the SU(2) isospin symmetric point, the spin 1/2 A resonance keeps having a mass close to the AJ/^p 
threshold, and ends up with a final mass of 4306.2 MeV (the AJ/ip threshold is at 4306.9 MeV). However, the spin 
5/2 and the two spin 3/2 states, that are also originated from this (164)2,0 C 662,0 pole, essentially disappear. One 
of the J = 3/2 states still shows up as a cusp very close to the AJ/ip threshold, and it has been included in the table. 
The second state with spin 3/2 (light magenta triangles in Fig. 3) and the spin 5/2 one appear as small unnoticeable 
peaks right at the Arjc and A J/\l/ thresholds, respectively. 

From the discussion above, the (164)2.0 C 662.0 HQSS multiplet could be incomplete. 

However the three A states (dark magenta circles for the two J^ = l/2~ states and dark magenta triangles for the 
J^ = 3/2~ resonance in Fig. 3) that stem from the (162)2,0 C 762, configuration of the light degrees of freedom 
turn out to be quite bound. Indeed, we find binding energies of at least 250 (150) MeV in the spin 1/2 (3/2) sector. 
These three states, nearly degenerate, form a clear HQSS multiplet. 

The models based on vector meson exchange, naturally predict a suppression factor in the baryon-meson amplitudes 
involving exchange of charm, from the propagator of the exchanged heavy vector meson. In the heavy quark limit, 
the suppression factor is of the order of l/mn-^'^ Therefore, in that limit, one expects a quenching of order Aly /Mo* 
for the charm exchanging amplitudes of those models. (Of course, the true factor for large but finite physical heavy 



The boson propagator is approximately 1/(2Mh{Eh — Mh)), with Mh the mass of the heavy vector meson and Eh its energy, and 
Eh — Mh is 0(1) in the heavy quark limit. 
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SU(6) X HQSS 
irrep 


SU(3) X HQSS 
irrep 


Mr 


Tr 


Couplings 
to main channels 


J 


702,0 


(102)2,0 


4005.8 


0. 


3Aj/^ = 0.3, S'E^B = 2.7, 5e,d' = 4.4, te'S* = 1-2 


1/2 


702,0 


(102)2,0 


4032.5 


0. 


Satjc = 0.2, QAj/i, = 0.1, gsjs = 2.9, gs.S* = 1-8 
<7s;£.* = 4.1 


3/2 


702,0 


(102)2,0 


4050.0 


0. 


ffA./zv. = 0.2, D = 0.8, 3e,d' = 1-9, 5e;d* = 5.1 


1/2 


562,0 


(104)2,0 


4306.2 


0. 

(cusp) 


3AJ/V. = 1-3, ffs.S = 0.3, 3e,B' = 0.3, 5e;d' = 0.3 


1/2 


562,0 


(104)2,0 


4306.8 


0. 

(cusp) 


ffA,,e ~ 0.1, 5a.7/^ = 0.8, gs^D = 0.2, = 0.2, 
SejS* = 0.1 


3/2 



TABLE IV: As in Table II, for the A ( J = 1/2, J = 3/2) resonances with hidden charm content. 



hadron masses needs not exactly coincide with this heavy quark limit estimate.) Our model is not directly based on 
exchange of vector mesons. Nevertheless, as commented above, we have verified that adding such suppression by hand 
in the charm exchanging amplitudes does not have an impact on our results. Even a factor (My/Mu*)^, proposed in 
the literature [68] has a very small effect in the position of the resonances we find. Presumably, this is due to the fact 
that the relevant channels have a small coupling. An exception comes from the two very weakly bound A resonances 
from the 562,0 ii'rep, which disappear due to the suppression of their dominant channel AJ/^p. 

V. CONCLUSIONS 

In the present work we develop a model for the interaction of lowest-lying 1/2+ and 3/2+ baryons with 0~ and 1~ 
mesons, including light and heavy flavors. The interaction is of zero range and it is modeled as a suitable extension of 
the Weinberg- Tomozawa term to include, besides chiral symmetry, spin-flavor symmetry in the light sector and heavy 
quark spin symmetry. These symmetries are only broken in our model to the extent that they are broken at level 
of the physical masses and meson decay constants. The OZI rule is also automatically implemented. Our extended 
WT model, Eq. (2.22), contains no adjustable parameters, although some ambiguity is present through the choice of 
renormalization prescription, as in all other hadronic molecular approaches. The model has been applied previously 
to the light sector and to charm or bottom sectors with a single heavy quark. Here we show that it admits a natural 
realization in sectors with hidden charm in such a way that HQSS is preserved. In particular, the spin of c quarks 
and the spin of c antiquarks are separately conserved. 

We have carried out a detailed analysis of the hidden charm sectors (i.e., with cc pairs) with C = 0,1,2,3 and 
their breaking as the symmetry is hfted from SU(6) x HQSS to SU(3) x HQSS (and then to SU(2) x HQSS and 
SU(2) of isospin). This allows to count the expected number of bound states or resonances, and to classify them into 
multiplets corresponding to the various symmetries. Taking the eigenvalues A's as undetermined (free) parameters, 
the results of Eqs. (2.30), (2.35), (2.38) and (2.42), for the C = 0, C = 1, C = 2 and C = 3 hidden charm sectors 
respectively, are general. Indeed, these equations fixes the most general structure of eigenvalues that can be deduced 
from SU(3) x HQSS. The rest of undetermined parameters, not fixed by this latter symmetry, accounts for non- 
diagonal transitions between multiplets with the same SU(3) x SUj^ (2) SF quantum numbers for the light degrees of 
freedom. Further, we have translated this general discussion of the group structure allowed by SU(3) x HQSS into 
Lagrangian form, for the charmless hidden-charm sector. This makes the HQSS of the model explicit and it allows to 
compare with other models in the literature. Finally, we have found the couplings of the HQSS effective Lagrangians 
of Eqs. (2.57)-(2.66) for the particular case of our extended WT model. This constitutes an additional check of its 
compatibility with HQSS. 

We have analyzed the charmless and strangeless sector, where we have dynamically generated several TV and A 
states. Actually, we predict the existence of seven A^-like and five A-like states with masses around 4 GeV, most of 
them as bound states. These states form heavy-quark spin multiplets, which are almost degenerate in mass. The N 
states form two HQSS multiplets. The lowest one has the light quark fiavor-spin content coupled to 82. Since the cc 
pair can couple to spin Sec = 0, 1, this HQSS multiplet consists of three nucleon states with J = 1/2, 1/2, and 3/2, 
and masses around 3930 MeV. On the other hand, the highest HQSS nucleon-like multiplet contains four resonances 
with J — 1/2, 3/2, 3/2 and 5/2, and masses around 4000 McV. In this case, these states are originated from the 
84 SF light configuration. These two SU(3) x HQSS multiplets arise from the 70-plet of SU(6) x HQSS. There are 
no N physical states coming from the 56-plet. With regards to A states, we find two multiples with very different 
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average masses, because in this case they are originated from different SU(6) x HQSS representations. The A multiplet 
coming from the (102)2,0 C 702,o irrep is formed by 3 states (J = 1/2, 1/2, 3/2) with an average mass of 4035MeV. 
Besides, we find only two (J= 1/2, 3/2) A resonances at the physical point out of the four states originated from 
the (104)2,0 C 562,0 hi the SU(6) limit. These two states are nearly degenerate, with a mass of 4306 MeV. 

When we compare our approach with the other molecular-type ones, we observe in our model a rich structure 
of resonances due to the many channels cooperating to create them. Our states are much more lighter that those 
predicted in the hidden-gauge scheme [74-76] , though significantly less bound that the crypto-exotic baryons reported 
in the zero range vector meson exchange model of Refs. [65, 66]. Moreover, we have presented the first prediction for 
exotic hidden-charm A-like resonances within a molecular baryon-meson scheme. 

In comparison with the quark model of Ref. [119], we find that our results are closer to those predicted by the FS 
hyperfine interaction discussed in Ref. [119], specially for the lowest lying states. 

The predicted new resonances definitely cannot be accommodated by quark models with three constituent quarks 
and they might be looked for in the forthcoming PANDA experiment at the future FAIR facility. 
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Appendix A: Spin-flavor states 

In this appendix we give details regarding the construction of the tensors AI'^b and B^^'^\ the SU(3) multiplets 
E°f,, Z)°, etc, and the computation of the matrix elements of the interaction. 

The wave-functions in spin-flavor space of the basic mesons and baryons are constructed in terms of bosonic quark 
and antiquark operators with spin and flavor labels, namely, Q'j^,Q^j^, ^ ^ u, d,s,c. The concrete wave- 

functions are those given in the Appendix A of [81] with the following modiflcation: a minus sign is to be applied to 
all 1/2+ baryons, to all 0~ mesons except 7?, r}' and ?7c, and to w and J/^ (denoted ip in [81]). No change of sign 
is to be applied to 3/2+ baryons, nor to ?7, 77' and ?7c, nor to 1~ mesons (except 4>, ui and J/ip)- 

The states just defined are standard with respect to the flavor and spin-flavor groups conventions of [111]. In 
particular they are SU(2)j, SU(2)/ standard and follow the convention of [116] for flavor SU(3) and SU(4). The only 
exceptions come from the neutral mesons for which we use ideal mixing. In terms of these, the standard states of 
[111] are given by: 

1 J?') Stan 
I??') Stan 

l^c)stan 

l^s) 

1^1) 

In these formulas, the right-hand sides contain the physical (or rather ideal mixing) neutral mesons that we use in this 
work. Their wave- functions are constructed as indicated above (i.e., from those in [81]). The left-hand sides contain 



= W) (SU(3)), (Al) 

= \IIW) + I\ric) (SU(4)), (A2) 

= -^l'/) + \/fl^c) (SU(4)), (A3) 

= + (SU(3) and SU(4)), (A4) 

= (SU(3)), (A5) 

= ^Jl\^)-l\'l^) + \\m (SU(4)), (A6) 
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the standard or mathematical states used in [111]. They have good quantum numbers with respeet to SU(6) or SU(8) 
(and their corresponding chain of subgroups). 

In order to construct the tensors M'^g and B^^'~^ , M^^b and B^^^j, with good transformations properties, the 
foUowing procedure is used: For aU flavors / = u, d,s,c, and for the various creation operators QJ^^, ^/i' ^^^i^ ^'^'^ 
Q^T. , appearing in the wave- functions of the hadrons, the following replacements are to be applied: 



Q}^^+Q}i, g}^^+Q}2, gt^^+Qt^ (A8) 

Note i) the minus sign in Q^j^, and ii) for quarks, the labels 1 and 2 correspond to spin up and down, respectively, 
but for antiquarks they correspond to spin down and up, respectively. 

After the replacement, there are only operators Q\, and Q^^ for creation (and Q^^ and Qa for annihilation), 
carrying any of the eight labels A = ul, dl, si, cl, u2, d2, s2, c2. These operators transform under SU(8) in the way 
indicated in Eq. (2.17). 

The meson matrix is then obtained by replacing Q^^Q\ with M^^ b and expressing it in terms of meson operators 
by inverting the wave-function equations. Similarly, for the baryons, Q\Q^bQ% '^^ replaced with -B^^p and then 
expressed in terms of baryon operators. The fields $"^3(3:;) and B^^'--^ [x) of Section JIB are constructed in the usual 
way from these annihilation and creation operators. 

The SU(3) multiplets introduced in Section II F arc obtained as follows. 

^abc _ j^{asi,bs2,csa)3/2 (A9) 

Here A — asi, B = 6.S2, C — cs^, with a, 6,c e {u,d, s}, si, 52,53 G {1,2}, and the notation {ABC)3/2 indicates that 
the spin part is coupled to J = 3/2 with 1 =t and 2 =|. Besides we refer here to the (annihilation) operator; the 
field AJ^''^(x) is constructed out of it. 



(AlO) 



Y^ab _ ^((asi,f)S2)icS3)i/2^ ^ (All) 
j^*a6 _ ^((asi,6s2)icS3)3/2 (A12) 



lg^j^^^((''Sl,CS2)oCS3)i/2^ • 



(A13) 



For the mesons 





— J" CS2)o' 




CS2) 


Vc 




J/4' 


— CS2) 



(A14) 
(A15) 
(A16) 
(A17) 

Recalling that Qfi = Qf^ and 0/2 -Qff. follows that Mt""^^^)^ = {M^'^bi + M°^62)/V2, while M^^^^,,)^ 
equals -M^^hz, {M^^bi - ^"^^2)!^ and iW^^^i, for J3 +1,0,-1, respectively. 

Finally, we remark that we systematically take the coupling of baryon and meson in the order | baryon) ® | meson), 
rather than | meson) ® | baryon). 

Appendix B: Baryon-meson matrix elements 

The coefficients D^j, appearing in Eq. (3.1), for the charmless (C — 0) and strangeless {S = 0) sector are compiled 
in this Appendix (Tables V, VI, VII, VIII, IX and X). In addition, we also provide here the corresponding coefficients 



20 Note that IV") of [81] corresponds to -| J/i/'} here, not to \^) of [111] and of Eq. (A7). 
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for the C = and 7^ sectors (Tables XI, XII, XIII, XIV, XV, XVI XVII, XVIII, XIX, XX, XXI, XXII). 
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TABLE V: C = 0, S = 0, / = 1/2, J = 1/2. 
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TABLE VL C = 0, S = 0, / = 1/2, J = 3/2. 
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TABLE VIIL C = 0, S = 0, / = 3/2, J = 1/2. 
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V 9 3 



TABLE IX: C = 0, S = 0, / = 3/2, J = 3/2. 
Arjc Alp E*5 EcL>* E*_D* 



A77C 
At/) 

EcL>' 
E!73* 






-2 







-V5 



2 



20 
3 

4 
3 

lU 

3 

20 
9 



TABLE X: C = 0, S = 0, J = 3/2, J = 5/2. 





Ai/) E*5* 




yT2 




VT2 4 



TABLE VII: C = 0, S = 0, / = 1/2, J = 5/2. 





E*D* 


E*£>* 


-2 



29 



TABLE XI: C = 0, S* = -1, / = 0, J = 1/2 



A?7c 
Aijj 
A^Ds 

AcD: 

EcD* 
E'^D' 
E'D' 






1 

1 

2 

'V3 






-^/3 



1 

V3 

^/2 





-^/2 
2 



\/2 

-1 





-1 

^/2 



-1 





-^/2 
^/2 



-\/2 V2 



6 

^/2 

-1 



2 
V2 



TABLE XIL C = 0, S' = -1, / = 0, J = 3/2. 

Atp AcD* EcD* EID E'^D* E*D* 



V2 sjl 
%/2 .[% 



-1 -1 - 



AV; 

Ac 5; 

E'D" 




2 
^/2 

2 
3 

10 
3 





-%/2 



= 0, s 


= -1, 








-2 





TABLE XVI: C = 0, S = -1, / = 1, J = 5/2. 





Jl*Tp E*I)* 






2 


^/8 




2 


^/8 




\/8 ^/8 


2 





TABLE XIX: C = 0, 5* = -2, / = 1/2, J = 5/2. 














VS 2 




V8 


2 x/8 




2 


Vs 



TABLE XX: C = 0, 5 = -3, / = 0, J = 1/2. 

nv- ^ -^1 

^'■ciJs \j 'J, y 3 3 y g 

o*n* -,/J -,fl ,fi _1 
"c-t^s y 3 y 3 y 9 3 



TABLE XXI: C = 0, S = -3, / = 0, J = 3/2. 
















V3 







%/3 


2 


-%/5 














-2 


/T 
V 3 


[20 

V 3 




-75 


[20 

V 


fi 
V 3 



-2 

4 
3 



10 
3 



TABLE XXII: C = 0, S' = -3, / = 0, J = 5/2. 





nip r2*_D* 




\/T2 




^/T2 4 



